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ABSTRACT 

Here we study the metallicity bias in the radial and tangential velocity dispersions, the de- 
rived quantity called anisotropy and the mean azimuthal velocity profiles of the MiUcy Way 
stellar halo using Blue Horizontal Branch (BHB) stars taken from SDSS/SEGUE survey. The 
comparatively metal-rich sample ([Fe/H] > —2) has prograde motion and is found to have 
an offset of 40 km s^^ in the mean azimuthal velocity with respect to a metal-poor sample 
([Fe/H] ^ —2) which has retrograde motion. The difference in rotation between the most 
metal-poor and most metal-rich population was found to be around 65 km s^^. For galac- 
tocentric distances r < 16 kpc, an offset in velocity dispersion profiles and anisotropy can 
also be seen. In the inner regions, the metal-poor population is in average tangential; however, 
anisotropy is found to decrease monotonically with radius independent of metallicity. Beyond 
r ~ 16 kpc, both the metal-rich and the metal-poor samples are found to have tangential 
motion. The metallicity bias in the kinematics of the halo stars qualitatively supports the co- 
existence of at least two-components in the halo having different formation history e.g. in-situ 
formation and formation by accretion. 
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1 INTRODUCTION 

The stellar halo is a n excellent resource for testing t heories 
of ga laxy formation jpreeman & Bland- HawthomI I2OO2I iHelmil 
|2008|) . With the advent of new high quality surveys involv- 
ing hundreds of thousands of stars, understanding the forma- 
tion of the halo has become even more challenging. Under the 
ACDM paradigm of hierarchical formation, the stellar halo is 
thought to have been produced in part by accretion of satellite 
galaxies. Observational evidence for this comes from the detec- 
tion of tidal streams and su bstructures in the halo of both the 
Milky Way and M3 1 (e.g., iMaiewskil 1 19931 : llbata et alj 1 19951 : 
iBelokurov et al Quantitative comparison of the amount of 

substructure and its physical properties seen in simulations with 
that of observations reveals broad ag reement w i th the accretion the- 
ory feell et al. 2008; Starkenburg et al.ll2009t [Gilbert et al.ll2009l : 
IXue et al]i201ll:ISharma et al.ll201 lb|). Recent cosmological simu- 
lations (e.g!. IZolotoy et al JI2009I : iFont et al.ll20 1 ll : iMcCarthy etsll 
l2012l : lTisseraetal.ll2012h including star formation and feedback 
suggest that the inner regions of the halo might be dominated by 
in-situ component which have chemical and kinematic properties 
different from that of the accreted component. The kinematic and 
spatial distributions of the in-situ component is thought to resemble 
that of the disk-flattening and net prograde rotation. 

Recent observational results also seem to suggest the need for 
an in-situ component. In their seminal paper, ICarollo et al. 1 j200l 
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hereafter C07) claim that the inner-halo component (r < 15 kpc) 
is comparatively metal-rich, has a prograde motion (0-50 km s~^), 
with a slightly flattened spatial distribution and with stars on mostly 
radial orbits. In contrast, the outer-halo is comparatively metal- 
poor, has a strong retrograde motion (40-70 km s^^), has a nearly 
spherical spatial distribution, with stars on a wide range of orbital 
eccentricities. These results were derived using stars in the Solar 
neighborhood which were roughly divided into inner and outer 
regions based on the maximum distance from the galactic plane 
reached by individual stellar orbits. Their analysis requires the use 
of proper motions for which it is important to have accurate dis- 
tances. The impact of dist ance errors on the retr ogra de signature 
has been debated, e.g. see lSchonrich et al] ( 1201 ih and lBeers et al.l 
(^01^. A cleaner distinction between the inner and outer regions 
can only be done by taking a sample of s tars which has much 
wider spatial distribution. This was done by iDeason et al.l ( 1201 lah 
using BHB stars grouped into four regions in metallicity and dis- 
tan ce spaces. Their resu lts are in qualitative agreement with CO? 
an d lCaroUo et al] ( l2010h . i.e., the metal-rich population is prograde 
while the metal-poor population is retrograde. However, they find 
that stars between 15-25 kpc are on circular orbits while those 
beyond 25 kpc are on radial orbits. The derived quantity called 
anisotropy dBinnev & Tremain3l2008h defined as 



/3 = 1- 



2(t2 



(1) 



describes the nature of the orbits, where ct^, o"9, and ct^ are radial, 
polar and azimuthal velocity dispersions respectively. Negative and 
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positive values of the anisotropy corresponds to the tangential and 
radial dominance of the orbits respectively. The arbitrary distance 
ranges make it difficult to compare the /3 results of BHB stars with 
those of C07. An arbitrary cut-off in metallicit y ([Fe/H] = —2) 
was used to group the stars. In our recent work dKafle et al]|2012l 
hereafter K12), we study /? as a function of radius in greater detail 
and we show that the profile is not a simple monotonic function. 
We find that for r < 15 kpc, stars are on radial orbits; at larger 
radii, there is a sharp drop in /3 with a minimum of /? « —1.2 at 17 
kpc; beyond 20 kpc, /? rises slightly. Given the complex variation 
of P with radius, it is important to study the radial profile of /3 and 
Vrot for different metallicities. 

This paper is organized as follows. In Section 2, we discuss the 
theoretical aspect of our analysis, and the data being used. Section 
3 presents our result about the metallicity biases in the kinemati- 
cal profiles, (Jr(r), CTe(r), o-0(r), /3(r), Vrot{r), of the Milky Way 
stellar halo. In Section 4, we conclude our result and describe the 
implications of the biases for the formation of the halo. 



2 THEORY AND ANALYSIS 

2.1 Distribution Function and Parameter estimation 

Our main aim is to estimate the rotational velocity and the velocity 
dispersion as a function of radius. To do this we model the veloc- 
ities using the Gaussian Velocity Ellipsoidal distribution function 
(GVE DF). A GVE DF with rotation about z-axis is given by 
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Figure 1. [Fe/H] distribution: Top-inset shows the error in WBG [Fe/H], 
quoted in the SEGUE database, and the bottom-inset shows the fraction of 
metal-poor stars with Poisson errors, both as a function of r. In the main 
plot, the black solid Hne is the [Fe/H] distribution of our sample, whereas 
the red solid line is the [Fe/H] distribution of the simulated stellar halo. 
The red dashed hne is the convolved red line (for [Fe/H]> —3.5) with the 
systematic error of 0.30 and the gaussian random eiTor of 0.40. 



as follows 



(4) 



p(r) 



■ exp 



(2) 



where v^ot is the mean azimuthal velocity. A GVE DF has been 
used in the context of the stellar halo to estimate the rotation 



I Frenk & White! I980f) and the velocity disp ersions dSirko et al] 
l2004bl ; ISmith et al.ll2009bl; iKafle et alj|2012l) . In their studies of 



halo subdwarf stars, ISmith et al. I (l2009ah and ISond et al.l ilOldt 



have found that the tilt of the velocity ellipsoid is small and consis- 
tent with zero, and is therefore ignored in this analysis. 

Ideally, we would like to know the full six-dimensional phase- 
space information of an individual halo star. But presently the 
proper motion information of the stars in the stellar halo beyond the 
Solar neighborhood (r > 10 kpc) is not accurate enough to prop- 
erly constrain the tangential motions. We thus resort to the line-of- 
sight component of the velocity and marginalize the DF over the 
two unknown components of the velocity (vi, vi,)- The marginal- 
ized DF can be expressed as. 



F(l,b, d, «ios \o-r,ae,a4,, v^ot ) 



f{r,\r)dvidvb. 



(3) 



Note the marginalization is over the tangential components of the 
velocities (vi, Vb) in the heliocentric rest frame, whereas DF (Equa- 
tion[2} is expressed in the galactocentric frame of reference. For the 
transformation between these reference frames, we assume the lAU 
adopted values of the Sun-Galactic centre to be 8.5 kpc and the cir- 
cular velocity at the Sun's position (ulsr) to be 220 km s^^. We 
also discuss the effect of changing the value of ulsr- The values 
for the peculiar motion of the Sun (U, V, W)^ = (+11.1, -1-12.24, 
-1-7.25) in km s^^ are adopted from lSchonrich et al.l(l2010l) . 

Using Equation [3] we can express the log-likelihood function 



where n is the total number of stars. Next we use the IVIarkov Chain 
Monte Carlo (MCMC) to compute the posterior distribution of the 
model parameters. The median values of the model parameters are 
quoted as our estimates, whereas 16*'' and 84*'' percentiles val- 
ues are used to compute the associated uncertainties. Note, in the 
GVE DF (Equation|2ll. the density term p{r) is not a function of the 
model parameters and thus has no effect on the likelihood analysis. 



2.2 Sample 

BHB stars are luminous and their distances can be measured with 
relatively more accuracy compared to other stellar populations, 
hence we use them to study th e stellar halo. Th e sample of BHB 
stars that we use is taken from lXue et al.l ( 1201 iL XII) comprising 
4985 BHB stars obtained from SDSS/SEGUE. To minimize con- 
tamination from disk stars, we select stars that are well removed 
from the plane of the disk (l^l > 4 kpc). We do not impose any 
kinematic limits to obtain our sample and hence for the purpose of 
kinematic studies this sample of BHB stars can be considered to be 
unbiased. To improve the distances, we recalibrate X II dist ances 
using a color-magnitude relation from Deason et al. ( 201 Ibt) . The 
dispersion in (?-band magnitudes of the data results a distance un- 
certainty of 6 per cent. In the radial velocity measurements, 95 per 
cent of our sample has an uncertainty of less than 9 km s^^. The 
final sample size of the BHB stars we use is 4386. 

The XI 1 catalog of BHB stars originally did not provide 
the corresponding metallicity values. To obtain [Fe/H], we thus 
cross-check SEGUE database for each star against the XII sam- 
ple. The [Fe/H] met allicity we adopt are the estimates based on 
IWilhelm et al. I( ll999h which is provided in the SEGUE database 
under the heading 'feh_WBG'. In Figure Q] we present the [Fe/H] 
distribution of our sample. In the inset of the figure, the errors in 
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WBG [Fe/H] estimates are shown as a function of distance (r) 
for botli the metal-rich (red line) and metal-poor (green line) sub- 
sample as well as for the total sample (black line). Note, errors in 
WBG [Fe/H] beyond r = 35 kpc are significantly larger (> 0.27 
dex). Also shown in the figure as a red l ine is the metallicity of BHB 
stars in the simulated stellar halo s by iBullock & Johnstonl ( |2005|) 
sampled using the code Galaxia i lSharmaetal.ll2011ah . The code 
Galaxia takes the N-body particles and their given age and metallic- 
ity distributions and then uses Padova stellar isochrones to generate 
stars. A sharp upper cut-off in the metallicity at around [Fe/H] ~ 
-1.7 can be seen and this constraint comes purely from the theory 
of stellar evolution. According to the synthetic isochrones used by 
Galaxia (Padova isochrones), there are no BHB stars with age less 
than 13 Gyr and metallicity greater than -1.7. In reality, globular 
clusters and the field stars (Behr 2003i) are found to have metal-rich 
BHB stars. We find that when the simulated BHB stars are con- 
volved with an uncertainty of 0.40 dex and an offset of 0.30 dex 
we can match the metallicity of the SEGUE BHB stars (dashed red 
line). The disagreement among the theoretical prediction and the 
observed metallicities of the BHB stars could possibly be either 
because the synthetic isochrone do not correctly model the BHB 
stars or may be because of inaccuracies in the computation of WBG 
BHB metallicities quoted in SEGUE database. 

Also shown in Figure [T] is the fraction of metal-poor stars 
([Fe/H] ^ —2) as a function of radius. It can be seen that the 
fraction increases with distance suggesting that the inner regions 
are slightly biased towards being metal enriched. However, in the 
outer parts, it starts to fall again due to the sharp increase in un- 
certainty of [Fe/H] at larger distances, e.g., the outer part could be 
metal-poor but a large uncertainty in [Fe/H] will spuriously classify 
metal-poor stars as relatively metal enriched. 



3 RESULTS: METALLICITY BIASES IN KINEMATICS 

We present here our estimates of the model parameters 
ar{r),ae{r),(T^{r), /?(r) and Vrotir) in comparatively metal-rich 
([Fe/H] > -2) and metal-poor ([Fe/H]^ -2) bins. To compute the 
kinematic profiles as a function of radius we choose a moving cen- 
tre scheme for binning. Briefly, in this binning scheme a set of equi- 
spaced positions in r are chosen and then at each position an equal 
number of points (jibin) from the central value are used to estimate 
the velocity dispersions. The number density of stars in r is not 
uniform and hence we report the mean r of the points in each bin 
as its final position. It should be noted that the bins are overlapping 
in this scheme; a detail discussion on associated pros and cons of 
using the moving centre bins are given in K12. 



3.1 Bias in Mean Azimuthal Velocity Profile 

We first compute the profile of mean azimuthal velocity (vrot in 
Equation I2J with nbin = 600. The blue markers in Figure |2] show 
the Vrot profiles of the total sample of BHB stars. It can be seen that 
the v-cot ir) profile of the overall population is nearly constant. The 
bottom panel shows the profiles for the metal-rich and metal-poor 
sample separately. In the inner region, r < 18 kpc, the metal-rich 
population is found to have a prograde motion of Wrot = 10 km s^^ 
which declines slightly in the outer-region and attains null rotation 
at r > 20 kpc. In contrast, the metal-poor population has a sub- 
stantial retrograde motion of 30 km s^^ for r < 18 kpc. At larger 
r, Wrot values of the metal-poor population also decrease to about 
18 km s"\ Fi gure|2]is restricted to r < 32 kpc as beyond here 
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Figure 2. Metallicity bias in Viot- 'Uiot estimates in two metallicity bins 
[Fe/H] > -2 (red) and [Fe/H] C -2 (green). The v^ot profile for the total 
sample is shown with the blue mai'kers. Black vertical solid line demai'cates 
bins that have stars with r < 30 kpc. 
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Figure 3. Urot dependence on [Fe/H] and the effect of tiLSH,: The blue dots 
with the error bars show the v^ot as a function of [Fe/H] for the sample 
of BHB stars, with i)lsr = 220 km s~^, whereas red open circles ai'e the 
values for the sample restricted to r < 30 kpc. Black and magenta mai'kers 
with the error bai"s are the corresponding estimates for vlsr = 200, 240 km 
s^^ respectively. 



the large uncertainty in [Fe/H] (see Figure [T) washes out the dis- 
tinction between the metal-rich and the metal-poor population and 
makes the v^ot profiles to converge to a common value. Addition- 
ally, the error in v^ot also begins to increase with r. There are no 
estimates for metal-poor samples for r < 12 kpc which is due to 
the fewer metal-poor stars in our sample in that regime than the 
met al-rich stars, as expect ed. 

iDeason et alj j201 lah also derive Viot and find that in a large 
bin ranging 10 < r / kpc < 50, the metal-rich population has a net 
prograde motion of 15 km s^^ whereas the metal-poor population 
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has a net retrograde motion of 25 km s~ .In agreement with Figure 
|2l it can be said that their estimates of «rot are more lilcely the better 
representation of the inner halo within r < 18 kpc, but definitely 
not of the halo out to r = 50 kpc. Note that their sample of BHB 
stars are not selected based on Balmer line-strengths, and thus we 
anticipate significant contamination from the Blue Straggler popu- 
lation as well as from t he cooler and low g ravity candidates of the 
horizontal branch stars jSirko et al.ll2004al) such as RR Lyrae. The 
fact that the metal-rich and metal-poor profiles do not have a very 
strong dependence on radius means that we can neglect the radial 
variation of v^ot and study the v^ot as a function of metallicity; this 
is shown in Figure |3] A switch from retrograde to prograde mo- 
tion at [Fe/H] ~ -1.8 is notable in this figure, suggesting a separate 
metal-rich component. We checked that this transition also remains 
for the sub-sample of stars within r < 30 kpc, for which eiTors in 
[Fe/H] is < 0.27 dex. Given about 0.22 dex uncertainity in [Fe/H], 
we find that a two component theory, Vrot being a step function of 
[Fe/H], can explain the above profile, but a smooth transition to the 
low metallicity end cannot be ruled out. From the fact that the v^ot 
profile is not symmetric about [Fe/H] — -2, we can infer that the 
actual transition to the low metallicity component will be at [Fe/H] 
< -2. This makes sense in the context of the metallicity distribution 
of BHB stars for the simulated stellar halo that was shown in Figure 
[T] namely that the synthetic isochrones predict very few BHB stars 
with [Fe/H] > -2. 

The ulsr has a correlation with Vrot estimates. To investigate 
this we repeat the analysis for different values of wlsr and the re- 
sults are shown in Figure|3] The effect of changing wlsr is to shift 
the curves up or down by the corresponding amount but keeping the 
feature of the profile nearly intact. Hence if vlsr is ~250 km s^^ 
(for Rq = 8.5 kpc: [McMillan & Binnevll201Cl) . then the metal- 
poor population will have nearly zero rotation and the metal-rich 
population will have a strong prograde rotation. This is consistent 
with the picture that the metal-poor stars that dominate the outer 
regions of the halo is formed by accretion of multiple satellites 
without any preferred direction of rotation. On the other hand, the 
metal-rich stars that dominate the inner-halo could have their ori- 
gin related to the disc through the in-situ stars, or due to a massive 
satellite which happened to have the rotation in the same direction 
as disc. However, irrespective of the adopted value of ulsr, the 
difference in the rotational properties of metal-rich and metal-poor 
population remains intact and our results show that the difference 
in Vrot between the most metal-rich and most metal-poor is at least 
65 km s~^. 

3.2 Bias in Velocity Dispersion Profile 

Apart from Viat, do the metal-rich and metal-poor subsamples dif- 
fer in any other kinematic property? To answer this, we study the 
metallicity dependence of ar{r),ag{r),a,f,{r) and /?(r), and are 
shown in Figure |4] The dispersions are measured in radials bins 
with ribin = 600 stars for r < 18.9 kpc and Tibin ~ 1200 for 
r > 18.9 kpc. This demarcation is shown by the vertical dashed 
line. We have to use larger ribin in the outer parts because fewer 
stars do not constrain the ae and a^. This is because in the outer 
parts the radial velocity carries less information about the tangen- 
tial components. A drawback of Tibin being large is that the spatial 
resolution is compromised. Nevertheless, large bins are useful to 
study general trends. In K12 it was shown that the overall trend in 
dispersion and /? profile of BHB stars remains intact even as ribin 
is increased from 750 to 1200, except for the fact that the profiles 
with large nbin are smoothed versions of those with lower ribin. 
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Figure 4. Metallicity bias in the velocity dispersions: red and green lines 
are the estimated values of velocity dispersions and coiTesponding /3 and 
a in the compai'atively metal-rich ([Fe/H] > —2) and metal-poor ([Fe/H] 
C —2) bins. Shaded area denote the errors associated. 
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Figure 5. Metallicity bias for the sample in a range 5 < r/kpc< 20 : the 
1(T and 2(T confidence contours for the joint-probability distributions of the 
parameters in two metallicity bins [Fe/H] > —2 (red) and [Fe/H] —2 
(green) respectively. 



From the figure we can see that below r = 16 kpc the trends 
in all the five listed profiles are similar for both the meal-rich and 
metal-poor populations. Namely, with increase in r, the Or de- 
creases, the (jg and cr^ increase and the /? decreases. However, 
the profiles are offset from each other. Also, the total dispersion 
a (= + Cfl + o"^)) is found to have an offset (40 km s^^). 

Overall the metal-poor population has lower /3 than metal-rich. The 
maximum difference is at r = 16 kpc, with ag being very high and 
the P being very low for the metal-poor population. Going from 
r = 16 to r = 20 kpc, one finds that the metal-poor population 
has almost constant (Jr whereas for the metal-rich population the 
Gr drops. The /? profile for the metal-poor population remains rel- 
atively constant whereas for the metal-rich it keeps on decreasing 
with minimum value being about -2. 

Beyond r = 20 kpc, the metal-rich and metal-poor profiles 
slowly start to converge. This is probably related to the fact that 
the uncertainty in [Fe/H] increases sharply beyond r = 30 kpc. 
So any kinematic differences between the metal-poor and metal- 
rich are expected to be washed out if stars beyond r = 30 kpc 
are included. Note in our scheme the bins are overlapping, so a bin 
at a given radius can contain stars within a wide range in r. To 
show this more clearly, in the figure we plot a solid vertical line 
that demarcates bins that have stars with r < 30 kpc. It can be seen 
that as we move to the right of this line the two cr,- profiles start to 
converge. 

The fall at r = 16 kpc in Or of the metal-rich population 
roughly coincides with the break radius in the density distribution 
of the halo stars; this fall was also noticed by Deason et a l. (2011ai) 
in the aias profile. In a recent paper dPeason et al.ll20I3h . it was 
postulated that a break radius in the density profile of a halo could 
be due to accretion of a massive metal-rich accretion event and a 
signature of this would be the fall in ar of the metal-rich popu- 
lation. We do see this dip in the metal-rich BHB stars but they 
are found to be on tangential orbits while shell-like structures in 
the lBullock & Johnstori feoOSi) simulations are generally found for 
satellites on highly radial orbits. It remains to be seen if a satellite 
on more milder radial orbit can still give rise to a break in density. 

It is apparent in the Figures |2]and|4]that the analyses are done 
in smaller bin sizes and consequently, the errors for each bin are 
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Table 1. Estimated parameters in the comparatively metal-rich and metal- 
poor bins in the range 5 < rA:pc< 20 for two different cases namely, when 
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larger. Thus, here we also provide additional measures of our model 
parameters, a, P and v-rot , for the metal-rich and the metal-poor 
samples covering the larger radial range 5 < r/kpc< 20. In Figure 
[5] we show the joint probability distributions of the model parame- 
ters for these sub-samples, where for the metal-poor and the metal- 
rich samples the confidence contours (Icr and 2cr) are shown in 
green and red colors respectively. The estimated values of these pa- 
rameters are given in Table[T] The differences among the estimated 
values of the kinematic parameters a and Viot of the metal-rich 
and the metal-poor samples are found to be ~ 50 km s^^. Also 
given in Table [T] are the estimates of the parameters for the sam- 
ple masked to remove the contamination from Sagittarius stellar 
stream and the Virgo overdensity . The cuts chosen to mask these 
structures are purely geometrical and are same as given in Section 
3.4 of K 12 (for the further details about the cuts see the references 
provided there within). It can be seen in the Table that the exclu- 
sions of these two substructures has negligible effect on the final 
results. 



4 CONCLUSIONS 

In this paper, we derived the mean azimuthal velocity, the velocity 
dispersion and the anisotropy profiles of the Milky Way stellar halo, 
using BHB stars, as a function of distance r and metallicity [Fe/H]. 

The comparatively metal-rich ([Fe/H] > -2) and metal-poor 
([Fe/H] ^ -2) populations of BHB stars within r < 16 kpc have 
distinct kinematic properties with offsets of 40 km s^^ in both Wrot 
and o. Additionally, an offset is also found in (Jr, og, and /3. 
This suggests the presence of at least two sub-populations in the 
halo. When v^ot is studied as a function of [Fe/H], we find that «rot 
decreases monotonically and then saturates to a constant value for 
[Fe/H] > —1.8. The saturation again suggests the existence of a 
distinct metal-rich component. The difference in rotation between 
the most metal-poor and most metal-rich population was found to 
be around 65 km s^^. 

Overall the /3 profile for both the metal-rich and the metal- 
poor populations show similar trends with radius, i.e., /? decreases 
with radiu s and then saturates a t a negative value. This is in agree- 
ment with lDeason et alj ( 1201 lal) in the sense that /3 is more a func- 
tion of radius than metallicity, but differs slightly from GOV who 
suggest that the metal-rich halo is radial while the metal-poor has 
a wide range of values of /3. Our results suggest that below r < 15 
kpc, the metal-rich population has radial orbits and the metal-poor 
population is on average tangential. But in the range 15 < r/kpc 
< 25, both the metal-rich and metal-poor populations are predom- 
inantly tangential. Finally, regarding the spatial distribution of the 
metallicity we find that although the inner regions (r < 15 kpc) are 
slightly biased towards being metal-rich, in general the metal-rich 
and metal-poor populations have substantial spatial overlap. 
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There are two factors which make it difficult to compare the 
BHB results with the findings of C07. Firstly, BHB stars are in gen- 
eral metal-poor so they preferentially sample the metal-poor popu- 
lation of the halo. Secondly, there is significant uncertainty in esti- 
mated metallicity. The uncertainty in fact increases sharply beyond 
30 kpc. The above two factors will lead to significant overlap be- 
tween the metal-poor and metal-rich populations and will smooth 
out the differences between the two populations. We do find evi- 
dence of this in the outer parts (r > 25 kpc) where a convergence in 
(Tr and Urot is noticed for the distinct populations due to the larger 
measurement uncertainties. Thus the physical differences between 
the two populations could be even stronger, but a deeper stellar sur- 
vey would be required to confirm this. 

Note a contamination of the BHB stars with other populations 
still a possibility. We find that theoretical isochrones do not predict 
BHB stars with [Fe/H] > —1.7 and have ages less than 13 Gyr 
A distribution of synthetic BHB stars can only be matched if we 
convolve with an uncertainty of 0.40 dex and allow for an offset of 
0.30 dex with respect to the metallicity of the SEGUE BHB stars. 
However, so me probable field BHB stars (e.g. BD+49°2137 from 
lBehj ( l2003h ) have been found to be metal-rich, thus it could be true 
that the synthetic isochrone do not model BHB stars well. Within 
r =16-20 kpc, we find that the cr^ and /? of the metal-rich pop- 
ulation decrease sharply while that of the metal-poor population 
show a slight increase. The fall in a of the metal-rich population 
roughly coincides with the break r adius in the density distribution 
of the halo stars. As suggested by iDeason et al.l ( 1201 3l) this could 
be a signature of accretion of a massive metal-rich accretion event. 

To conclude, we find that the metal-rich and the metal-poor 
populations of the halo BHB stars have distinct differences in kine- 
matic properties. But note that there is no sharp transition from one 
population to another e.g. C07 [Fe/H] distributions of the inner- 
halo and the outer-halo peak at -1.6 and -2.2 respectively but have 
a spread which are larger than the difference of their peak values. 
Therefore a simple division into two groups by [Fe/H] will include 
substantial numbers of the outer-halo in the metal-rich group and 
vice-versa. 

The difference in the metal-rich and metal-poor seen by us is 
in qualitative agreem ent w ith earlier obs e rvation al findings of C07, 
ICarollo et al.l ( |2010|) . and iDeason et al.l ( 1201 lal) and supports the 
dual formation scenario as seen in cosmological simulations, i.e., 
the inner halo is dominated by an in-situ component and the o uter 
halo by an accreted compon ent (e.g.,iZolotov et al. 2009; Font e t al.l 
|201 iMMcCarthv et alj20l3) . However, to accept or reject the above 
hypothesis, a detailed comparison of observations with simulations 
still needs to be done. Also a comparative study of the kinematics 
of stellar populations other than BHB stars such as turn-off stars 
(fairer sample in metallicity) or giants (probing farther region of the 
outer-halo) or standard candles e.g. RR Lyrae (accurate distances) 
is needed. 
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